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Abstract 

We study and compare the decoherent histories approach, the environment- 
induced decoherence and the locahzation properties of the solutions to the 
stochastic Schrodinger equation in quantum jump simulation and quantum 
state diffusion approaches, for a quantum two-level system model. We show, 
in particular, that there is a close connection between the decoherent his- 
tories and the quantum jump simulation, complementing a connection with 
the quantum state diffusion approach noted earlier by Diosi, Gisin, Halliwell 
and Percival. In the case of the decoherent histories analysis, the degree of 

approximate decoherence is discussed in detail. 
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I. INTRODUCTION 



Two primary paradigms - the environment-induced decoherence approach, proposed by 
Zurek and the consistent histories approach by Griffiths and later by Omnes 

and by Gell-Mann and Hartle have been recently developed to solve the fundamental 
issues in quantum theory, especially, quantum measurement problems and the transition 
from quantum to classical. The environment-induced decoherence emphasizes the division 
between the system and its environment. The interaction of the system with its environment 
is responsible for the decay of the quantum coherence of the system. The decoherent histories 
approach is designed to provide the most general descriptions for a closed system by using 
the concept of history - a sequence of events at a succession of times. Both approaches are 
applicable to open quantum systems. 

Another set of viable theories within the framework of quantum mechanics are the various 
unravellings of master equation as stochastic Schrodinger equations for the single member of 
the ensemble. Among others, the quantum state diffusion and the quantum jump simulation 
approaches have been extensively studied in recent years {e.g., see P-[12|). As phenomeno- 
logical theories, these stochastic approaches are not only of theoretical interest but also of 
practical value. 

The open quantum system provides a unified framework to exhibit the properties of the 
various approaches we have mentioned above. The master equation which describes the evo- 
lution of the open quantum system plays a central role in the investigations into the decay of 
quantum coherence due to the interaction with a much larger environment. However, it does 
not tell us how an individual member of an ensemble evolves in a dissipative environment. 
The unravelling of master equation as stochastic Schrodinger equation could provide such a 
description within its domain of applicability. Corresponding to the decoherence process in 
the density operator formalism, in stochastic Schrodinger equation approaches, the solution 
to the stochastic Schrodinger equation often possesses a very remarkable property - the 
solution tends to localize at some special states after a localization time scale. For quantum 
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diffusion approach, this locahzation property has been justified in many different situations 

EMI- 

It is worth emphasizing that a key point in these approaches is the mutual influence 
between the system of interests and its environment. This mutual influence is the com- 
mon sources of the many different phenomena such as dissipation, fluctuation, decoherence, 
localization, etc. 

Analysis of decoherence and localization properties is usually rather involved. The en- 
tanglement of the complicated mathematics and the subtle conceptual issues often tends to 
make the detail scrutiny of the basic concepts impossible. The attractiveness of two-level 
system model is that, perhaps it is one of the simplest yet physical meaningful models. 

The purpose of this paper is to employ a widely used two-level system model as a unified 
framework to examine the dynamics of the open quantum system by the decoherent histories, 
the environment-induced decoherence and the stochastic Schrodinger equations. We mainly 
consider the decoherence process and the localization process. The various time scales 
concerning these processes are discussed. We provide a detailed analysis of the degree of 
decoherence which is important for a real physical process. 

The plan of this paper is as follows. In Section II, we briefly present our two-level system 
model and its basic properties. We study the consistent histories approach, the degree of 
decoherence in Sections III and IV, respectively. We study the unravelling of master equation 
and the localization properties of the solutions to the stochastic Schrodinger equations in 
both quantum jump simulations and quantum state diffusions in Section V. We make some 
remarks in Section VI. In Appendix we present a proof of Theorem in Section IV. 
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II. THE MODEL 



A fundamental building block in quantum theory is the two-level system. We consider a 
two-level atom system, which is radiatively damped by its interaction with the many modes 
of a radiation field in thermal equilibrium at temperature T. The upper level and lower level 
are denoted by |2) and |1), respectively. A widely used master equation for the two-level 
atom takes standard Lindblad form (in the Schrodinger picture) {e.g., see |TD|jn[]): 



P 



7_, 



+ ■^(n + l){2apa^ 



a^ap — pa) a) 



-f —n{2a^ pa — aa^ p — paa^) . 

Here, the Hamiltonian of the atom in the absence of the environment is given by 

H = cr,. 
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(2.1) 



(2.2) 



The Lindblad operators, which model the effects of the environment in this situation, are 

Li = J^(n + l)a, L2 = v^a^ (2.3) 



/^yna' . 

The transition rate from |2) |1) is described by the term proportional to {'~f/2){n + 1), 
and the transition rate from |1) |2) is described by the term proportional to (7/2)7?,. We 
use (Tx,(Jy and cr^ to denote Pauli matrices and a, a''" atomic lowering and raising operators, 
which are defined in the usual way 
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(2.4) 



and 
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2' 



= -{(^x- io-y), a^ = -{a^ + iay) 



(2.5) 



The master equation ( |2.1|) in the basis |2), |1) may be written as 
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Pii = -l{n + l)pii + 7np22 , 
P22 = j(n + l)pu - 7np22 , 

P12 = (2"- + ^) + MP12 , 

P2I = + ^) ~ ^^I'^^l • 



(2.6) 
(2.7) 
(2.8) 
(2.9) 



The general solutions to Eqs. (|2.6| )- (|2.9| ) are as follows 



Pii(t) = ^5i + 52e-^(^^+^)* 
n + 1 

P22(t) = 5i-i?2e-^(2^+^)S 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



where Bi{i = 1,2,3,4) are arbitrary constants which can be easily determined once the 
initial condition is given. For the initial density matrix with Tr(po) = 1 we easily get 

n + 1 



2n + l 



n 



B2 = pii(O) 

^3 = Pl2(0) 
54 = P2l(0) 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



In the basis |2), |1) the off-diagonal elements tend to zero as t goes to infinite. The 
density operator p tends to the stationary density operator p^: 
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2fi+l 



(2.18) 



This is an elementary example of environment-induced decoherence. Diagonalization occurs 
in the basis |2), |1). From equations ( |2.8| ) and ( |2.9D , it follows that the decoherence time 
scale to is given by 



It is easily shown that master equation is invariant under unitary transformations 
of the Lindblad operators: 

a I — > UaU\ a) i — > Ua)U\ (2.20) 

where ?7 is a unitary matrix. Correspondingly, the density operator p transforms in the 
same way: 

p I — > UpUK (2.21) 

Thus, when t — > cxo 

p — > UpsUl (2.22) 

Generally, the density matrix U psU'^ is no longer diagonal. This indicates that environment- 
induced decoherence does not occur in other bases. 



III. CONSISTENT HISTORIES IN THE TWO-LEVEL SYSTEM MODEL 

The decoherent histories approach [^H^ offers a sensible way to assign probabilities to a 
sequence of properties of a quantum system without referring to the measurements or to a 
classical domain. 

A history is defined in general as a sequence of properties of a closed system occurring 
at different times, which is denoted as 

c. = Pi?(U,-,^£Hti), (3.1) 

where Pj^^{ti) are the projection operators in the Heisenberg picture at times U: 

here, H is the Hamiltonian of the closed system. These projection operators satisfy exhaus- 
tive and exclusive conditions: 
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E^?(^0 = I. pSit^p^^t^ = s^^p^pjS{t^. (3.3) 

ai 



The superscript (i) labels the set of projections used at time ti and ai denotes the particular 
alternative. 

A natural way to assign the probability to a history is 

p{C^) = Tr(C,p(to)Ct) 

= Tr (Pi? (tn)...P«(ti)p(to)P«(ti)-^i?(tn)) • (3.4) 



However, one finds that ( |3.4| ) generally does not satisfy the usual probability sum rules. The 
necessary and sufficient condition to guarantee that the probability sum rules hold is that 
the real part of the decoherence functional D{a,a') vanishes, 

ReD(a,a') = ReTi{Cap{to)Cl) = (3.5) 



for any two different histories Ca and C^/. The sets of histories satisfying ( P3| ) are said 
to be consistent (or weakly decoherent) . Physical mechanisms causing ( pj.Sj ) to be satisfied 
typically lead also to the stronger condition 

D{a,a') = 0, Va ^ a' (3.6) 

which is called medium decoherence (In this paper, we simply refer it as decoherence) 

Although the decoherent histories approach was primary designed for a closed system, 
the approach is of particular importance for the open system which may be regarded as a 
subsystem of a large closed system. For the open quantum system, a natural coarse-graining 
is to focus only on the properties of the distinguished system whilst ignoring the environment. 
In this case, a natural selection of projections at each time is of form ® where 
is a projection onto the distinguished subsystem and denotes the identity projection on 
the environment. In the Markovian regime, the decoherence functional could be constructed 



entirely in terms of the reduced density matrix of the system [lIS 



= Tr iP!C^Kl_,,[Pt-l^...Kll[P!^^ , (3.7) 



where the trace is taken over the distinguished system only. The quantity KII_^ [ ' ] is the 
reduced density operator propagator: pt = Kq[pq]. 

In what follows we shall make a detailed analysis of the decoherent histories in the 
two-level model described by the master equation ( |2.1|) which depicts a Markovian process. 

First, let us consider the projection operators represented by 

Pi = |1)(1| and P2 = |2)(2| . (3.8) 

Obviously, {Pi, i = 1,2} form a set of complete and exclusive projection operators. Phys- 
ically, Pi may represent that the atom emits a photon whereas P2 may represent that the 
atom absorbs a photon. Then the decoherence functional at two time points is given by 



D[a,a:] = 5,,,^,Tr (P^^^ [P^K',^ [po] P^,\) . (3.9) 

It is easily shown that, for any 2x2 matrix A, the matrix PiAPj {i ^ j) is an upper (or 
a lower) triangle matrix. From equations ( p.8| ), (|2.9| ) we know that Kj'.^^l ■ ] propagates the 
matrix with zero diagonal elements into the matrix with zero diagonal elements. So, for 
any initial density matrix pq, the trace in Eq. ( |3.9| ) is exactly zero for any different pairs of 
histories {a 7^ a') and for any interval t2 — ti. This demonstrates that the set of histories 
consisting of projectors ( |3.8| ) are exactly decoherent. The generalization to n time points 
is straightforward. The exact decoherence for any time interval is slightly surprising. (The 
density matrix, by contrast, only becomes exactly diagonal as t ^ 00). This exactness is 
due to the simplicity of the model and we do not expect it to be a generic feature. 

Next, consider more general projection operators which correspond to the projection to 
any direction. With any direction denoted by a unit vector n = (sin 6 cos 0, sin 6 sin 0, cos 0), 
we associate a vector \n) which belongs to the Hilbert space of the two-level system, 

|n) =cos^|l) -e-^'^sin^|2) . (3.10) 

Then one can define the following projection operators on the Hilbert space of the system: 

P+=\n){n\, P. = \n){n\, (3.11) 



where \n') is the orthogonal comphmentary of \n): 

= e^-^sin^ll) +cos^|2) . (3.12) 

We shall show that a set of histories consisting of the projection operators P4. and P_ 
are approximately decoherent. To this end, first, note that for any 2x2 matrix A, 

Tr {P+AP_) = Tr (P-AP+) = . (3.13) 

Hence, from ( p.lO|) -( ^.13| ), it can be seen that, after the propag ation of Kil+'[-], all of the 
diagonal elements of matrix Kl'^[P±Kl^^[pQ]Pzf:] contain an exponential damping factor 

Damping factor = 6"^^^"+^)* . (3.14) 

Thus, we conclude 

~ 0, Va ^ a'. (3.15) 

This proves that the set of histories consisting of P+,P_ are approximately decoherent if 
time interval between tk and tk+i is larger than the characteristic time scale, 

^decoherence ~ /„ | 7T- (3.16) 

7(2n + 1) 

This is an expected result. We will give a more detailed estimate of the degree of decoherence 
in the next section. 

Note that tdecoherence dccrcascs as the coupling 7 is made stronger. To see the infiuence 
of the temperature of the bath, a simple expression for n may be chosen as 

From this, it is easy to see that the decoherence is more effective if the temperature of bath, 
T, increases. Conversely, decreasing temperature will make system spend more time to de- 
cohere. The maximum decoherence time scale is I/7 which corresponds to zero temperature 
of bath. 



In summary, we find exact decoherence of liistories cliaracterized by tlie projections onto 
|1) and |2), and approximate decofierence in any otlier basis. 

Finally, we examine tlie probabilities for two times histories consisting of projections 
(p.SI). These are given by 



P(l,2) 



P(2,l) 



P(2,2) 



[P,Kll[P,Kll[p,]P, 
n + 1 



;i-'5)pii(t 



1)^ 



2n + 1 

It (P,Kll[P,Kll[p,]P, 



n 



[l - S)p22{t\ 



2n + 1 

Tt(p^k\i\p^k\i\p,\p^ 
1 



2-^^ [^+(^+l)'^]Pii(^i), 
Tr(P2<nP2i^i*;NA]) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



Where 



Pii(^i' 



n 



2n + l 
71+ 1 
271+1 



Pii(O) 
Pii(O) 



n 



2n + l 

n 



(3.22) 

(3.23) 



271+1. 

and 6 = exp{— 7(271 + l)At} (At = ti — ti_i,i = 1,2). As for the n times histories, the 
calculation for the elementary probabilities will be straightforward. For instance. 



p(l, 1, ■ ■ ■ 1) = Tr {P,Kl:_^ [P.Klzl [Pi ■ ■ ■ Kll [P.Kll [p,]P, 
1 



Pi]Pi] 



= 7 [77 + (77 + 1)5]""' pn(ti) 

Similarly, one may calculate the transition probabilities, etc. 



(3.24) 



IV. DEGREE OF DECOHERENCE 



Physically, one would not expect the decoherence takes place exactly. Therefore the 
investigation of the approximate decoherence is of importance. In practical problems, one 
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can, at best, only expect that probability sum rules are satisfied up to order e, for some 
constant e < 1. Namely, the interference terms do not have to be exactly zero, but small 
than probabilities by a factor of e. One simple inequality which turns out to be very useful 
to the study of the degree of decoherence is pO| , |2T| : 

< t^D[a.(^D[d.d] ■ (4.1) 

We say that a system decoheres to order e if the decoherence functional satisfies ( ^111) . As 
shown in [^, such a condition implies that the most probability sum rules will then be 



satisfied to order e. 

Based on this two-level model we will study the degree of decoherence in some detail. 
To begin with, we establish the following trace inequality which is useful to our studies of 
the approximate decoherence. 

Theorem: Suppose that M and are two n x n positive definite matrices. Let P and Q 
be two n X n Hermitian matrices satisfying 

QP = PQ = 0. (4.2) 

Then 

|Tr {MPNQ)\^ < e^Tr (MPNP) Tr (MQNQ) (4.3) 

where e = min{e^, e^}, here = (A^, - X^J/iX'J^^ + X^J, = (A^,, - A^,J/(A^,, + 
■^min)) -^max) "^min -^max) -^min the maximal and the minimal eigenvalues of M and N, 
respectively. 

Remark: In fact, the condition that both M and are the positive definite matrices 
could be generalized to that one is positive definite, say M, while another A^ is positive 
semi definite. In this case, e = e*^. It is hoped that the above theorem is also useful in some 
other cases. The theorem is proved in the Appendix. 

For a general initial state represented by po (pure or mixed state), the decoherence 
functional of two time points may be written as 
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«'] = Tr (P±<^ [P.Kl^ [po] P- 



(4.4) 



We now write 



(4.5) 
(4.6) 



Then Eq. ( [4.4|) may be rewritten, in the new notation, as 



D[a,d] = Tr(5P_AP+). 



(4.7) 



Note that K in ( |4.6|) is the super-propagator for the projection operators 



(4.8) 



The evolution equation for the projection operators is given by 



7, 



+ -^{n + l){2a^Pa - a)aP - Pa^a) 
+ '^n{2aPa^ — aa^ P — aa^P), 



(4.9) 



where H, a, are defined as before. Note that the evolution equation for the projection 
operator P is different from that for the density operator p ( |2.1|) . This reflects the difference 
between the Schrodinger and Heisenberg pictures in the density operator formalism. The 
explicit form of Eq. ( |4.9| ) may be written 



Pii = 7(n + l)(P22-Pii), 

P22 = in{Pu - P22) , 

Pi2 = -[|(2n + l)-MA2, 

P21 = -[|(2n+l)+HP2i. 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



The general solutions to the above equations are 
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P22{t) = C,- -^C2e-^(^^+^)* , (4.15) 
n + 1 



Pi2(t) = C3e-[2(2"+i)-Hi ^ (4.16) 

P2l(t) = ^^e-[l(2n+l)+*a;]t ^ (4.17) 

where {i = 1, 2, 3, 4) are arbitrary constants. For given initial values, these constants can 
be expressed as 

C.^^Pn(0) + ^P.(0), (4.18) 

C^2 = ^^(Pii(0)-P22(0)), (4.19) 

Ca = Pi2(0) , (4.20) 

C4 = P2i(0). (4.21) 



From the definitions ( ^.5| ) and ( [4.6| ), it is easy to see that in general, both A and B could 



be positive definite matrices, and since P_ and P+ are projection operators, so the condition 
( [4.2| ) is automatically satisfied. Using the theorem above, we immediately arrive at 



|Tr {BP+AP_) |2 < e^Tr {BP+AP+) Tr {BP^APJ) . (4.22) 



That is, 



D[a,d]\ < ^ D[a,a])D[a!,a'] . (4.23) 



where e = min{e , e } 



e^ = \Xt-X^\, (4.24) 
= ■ (4.25) 



Here Xf (z = 1, 2) and Af (i = 1, 2) are two eigenvalues of A and P, respectively. (Note that 
+ = 1) From Eqs. ( [4.24| ) and ( [4.25| ), it is easily seen that the degree of decoherence 



may depend on both the projection operators we use and the initial state of the system. This 
is also an expected result. For the two-level system, and can be calculated exactly. 

13 



Consider, first, the eigenvalues of A. Since Eq. ( |2.1| ) preserves tlie trace, can be written 

as 



(4.26) 



Tlie determinant of A can be explicitly evaluated from the general solutions ( p.lO| , |2.11 
232| , |l3D, 



n 



rr+ 1 



2rr(rr+ 1^ 



+ 



Pii(O) 



2n + 
n 



Pl2(0)p2l(0) 



2n + 1 



P22(0) 



n + 1 



52 



6 



(4.27) 



2n + IJ 

where 5 = exp{— 7(2?!+ l)ti}. In order that a set of histories are decoherent, one expects 
that 6 should be small. 

Similarly, can be expressed as 



4det5 



'1 



V (Trfi)2 • 

From ( |4.14| )-( p:.17| ), TtB and deti? can be easily obtained: 



(4.28) 



TtB = 2Ci + ^C2S^, 
n + 1 



det B = Cl + 



— — -C1C2 — C3C4 
n + 1 



n 



5l - ^^2^?, 



(4.29) 
(4.30) 



where (5i = exp{— 7(2?!+ I)(t2 — ^i)}- The above discussions show explicitly how the degree 
of decoherence is related to the projection operators, the initial states and the temperature 
of bath, as well as the time spacing interval, in accordance with our general expectations. 

In the long time limit, the density matrix will tend to the stationary density matrix. 
Then we may get a much simpler expression for e^: 

1 



2n+ 1 



(4.31) 



As mentioned before, for the decoherent histories, S and 5i should be small. If we only keep 
the terms up to the first order of 61, then becomes 
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The expression for can be obtained from Eq. ( 4.27|) . 



1/2 

(4.32) 



It is seen from the above expressions that the degree of decoherence improves as the 
bath temperature increases. We also see that the projections with the smaller off-diagonal 
elements will give a better degree of decoherence. For a given system with the initial state, 
then the matter for investigation is to determine which histories, i.e., which string of pro- 
jections, will lead to the decoherence condition being satisfied. Therefore, we see that 
serves as the main criterion for the degree of decoherence. 

It is also of interest to compute the von Neumann entropy of pit) P,P2|. We find that the 
von Neumann entropy supplies a restriction on the degree of the decoherence of histories. 
It exhibits a tension between the predictability of the quantum state and the degree of 
decoherence. 

The von Neumann entropy provides a convenient measure of the loss of predictability: 

5* = -Tr(plnp) , (4.33) 

By definition, the more predictable state (pure state) may have less increase of the entropy 
in a fixed time period. This characterization process of predictability is called the "pre- 
dictability sieve" (coined by Zurek) which has been studied recently in quantum Brownian 
motion model by using the linear entropy P, P^P3| . We will see that two- level system pro- 
vides a very simple model to exhibit the relationship between the degree of decoherence and 
entropy by directly using the von Neumann entropy. 

For the purpose of the evaluating the entropy, we choose a special basis in which p is 
diagonal. Let Ai and A2 be the eigenvalues of p, then Eq. (|4.33|) reduces to 

S = -j2\\^Xi. (4.34) 



i=l 



Obviously, Ai and A2 can be expressed as 



1 + 1 - 

Ai = A2 = ^- . (4.35) 
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Hence, Eq. ([4.34|) can be rewritten as 



S 



[1 + 6^) (1 + 6-^) ^ (1 - e^) (1 - e^l 



(4.36) 



2 2 2 2 

It can be shown that >S'(e"^) is a monotonically decreasing function of e"^. Here, we find 
that there is a tension between the predictabihty and the degree of decoherence. That is, 
the initial density matrix which leads to less entropy production will give worse degree of 
decoherence, namely, would be not small (Here we are not considering e^). This tension 
between the predictability and the degree of the decoherence is a physically expected result. 
To obtain the higher degree of the decoherence one would expect that the environment has 
stronger influence on the system of interest, such as increasing the temperature of the bath. 
Then the predictabihty of the state, correspondingly, decreases. 

As an example, it is easy to see that the pure states |1) and |2) will lead to the largest 
entropy production. However, we have seen that the histories consisting of the projection 
onto these two levels give the best degree of decoherence (In this case is zero). In fact, it 
has been shown that the degree of decoherence is related to both the initial density matrix 
and projection operators used in the histories. 



V. UNRAVELLING OF MASTER EQUATION 

The master equation provides an ensemble description of a quantum system. The un- 
ravelling of master equation as the stochastic Schrodinger equation for the state vector has 
provided many insights into the foundation of quantum theory, especially in quantum mea- 
surement and the useful tools to study various practical problems in the quantum optics {e.g., 
l^,|l^]). In this section we will study the localization in the two different unravellings of 



see 



the master equation - quantum jump simulation and quantum state diffusion approaches. 
The former use the discrete random variables whereas the latter use the continuous random 
variables. 
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A. Quantum Jump Simulation 



The stochastic Schrodinger equation used in the quantum jump simulation describes 
a single quantum process. The system undergoes a smooth evolution until a jump takes 
place. This jump process is characterised by the real random variables which only take the 
discrete values. For the master equation (|2.lD, the stochastic Schrodinger equation takes the 
following form: 

\d%l)) = ~^H\iP)dt 

Here Li = \J'~^{ji + L2 = a/tWo^ are the Lindblad operators representing the influence 
of the environment and Ni = L\Li [i = 1,2). (Ni) = {iplNilip) represents quantum average 
and M represents the ensemble average. The real random variables dWi {i = 1, 2) satisfy 

dWidWj = SijdWi , (5.2) 
M{dWi) = {Ni)dt {i = 1,2). (5.3) 

Under condition (|5.2|) , it is easy to see that dWi only take two values: and 1. The master 



equation (|2.1| ) can be recovered from the stochastic Schrodinger equation ( |5.1|) in the sense 



that if l^p) is the solution to Eq. ( ^.1|) then p = M\ilj){ip\ satisfies master equation (|2.1|) . 



In what follows we shall discuss the the 'localization' properties of the single jump tra- 
jectories. Here, by 'localization' we mean that the quantum state vector generated by the 
stochastic Schrodinger equation will converge to some fixed states in the mean square . The 
physical meanings of this localization will become clear from the later discussions. 

More precisely, let A be an operator (not necessarily Hermitian), then we define the 
quantum mean square deviation as 

a{A,A) = {A^A) - {A^){A). (5.4) 
17 



If the solution of the stochastic Schrodinger equation (|0| ) satisfies 

M^a(AA)<0, (5.5) 

namely, the dispersion of the operator A tends to decrease as time evolves. Then we say 
that the solution localizes at the eigenstates of the operator A {A is sometimes called the 
collapse operator). 

For the stochastic Schrodinger equation for the quantum jump simulation in two-level 
system, the collapse operator is az- Then quantum mean square deviation in this case is 

(Aa,)2 = 1 - (a,)2 . (5.6) 

In order to prove the localization, we should first derive the evolution equation of the ex- 
pectation value of az by using the following formula: 

d{A) = {tlj\A\d^) + {dtlj\A\tlj) + , (5.7) 

where A is an operator. From (|5.1| ), it is straightforward to arrive at the following equation 

d{az) = (1 - {az))dW, - (1 + {az))dW2 

+ [{az) {m + N2) + {N, - N2)] dt . (5.8) 

Notice that 

diAazY = -2{az)d{az) - {d{az) f . (5.9) 

Then, inserting Eq. (|5.8|) into the above equation, taking the ensemble means and remem- 
bering ( ^.3|) , we obtain 

M|(Aa,)^ = -]p{n + l){l - {az)f{l + (a.)) 

-l^n{l + {az)f{l-{az)). (5.10) 

The right-hand side of Eq. ( |5.10| ) is non-positive, and that it vanishes if and only if is 
1 2) or |1). Hence we conclude that the solution to the stochastic Schrodinger equation ( ^.1|) 
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will localize at |2) or |1) after a certain time. That is, any initial state (which will be a 
superposition of |1) and |2)) will tend to a solution in which the atom undergoes stochastic 
jumps between |1) and |2). 

Let us now estimate this localization time. From ( p.lOp , a few manipulations directly 
give 

M^(Aa,)2 < -7(2n + l)(Aa,)2 . (5.11) 
So the localization rate tiocaiization is 

1 

^localization ~ /„ ; 7T ; (^•-'-^) 

7(2n + 1) 



which agrees with the decoherence time scale ( p.l9|) . 



In order to see the meanings of the localization of the quantum jump process, let us 
compute the evolution of the populations of the two levels: 

d\{im' = ^{{a,) - l)dW, + ^(1 + {a,))dW, 

-\ [(^.) {N, + iVa) + (iVi - N^)] dt , (5.13) 
rf|(2|^)P = i(l - {a,))dW, - \(\ + (a,))rfl^2 

+ ^ [(^.) (A^i + ^i) + (iVi - iV2)] dt . (5.14) 

Then it follows that from the above equations, 

IdWP^^. (5.15) 

That is to say, due to the influence of the environment to the system, after a localization time 
scale ( |5.12 ), the average populations of the first and second levels will become constant. This 



may assist to understand the meanings of the localization in the quantum jump simulation. 

In some sense, that the localization in quantum jump simulation chooses the basis |1), |2) 
appears to be natural, since they correspond to the trajectories that would actually observed 
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in an individual experiment. As expected, the set of histories consisting of projection onto 
the basis give the best degree of decoherence. In addition, we have seen that density matrix 
become diagonal in this basis. Here, we have demonstrated a close connection between the 
different approaches. 

B. Quantum state diffusion 

In this subsection, we will illustrate the localization process in another unravelling of the 
master equation - the quantum state diffusion approach, which was introduced by Gisin and 
Percival to describe the quantum open system by using a stochastic Schrodinger equation 
( which is often called the Langevin-Ito stochastic differential equation) for the normalized 
pure state vector of an individual system of the ensemble. Similar to the quantum jump 
simulation, a solution of the Langevin-Ito equation for the diffusion of a pure quantum state 
in state space represents a single member of an ensemble whose density operator satisfies 
the corresponding master equation. 

Generally, if the master equation takes the standard Lindblad form: 

p = -'-[H,p\+Y.{UpL\ - h\Up-]^pL\U) (5.17) 

Then, correspondingly, the Langevin-Ito stochastic equation can be written as 

\dip) = -lH\il))dt 
n 

+ T.i{Ll)L.-lLlL,-l{Ll){U)mdt 

i 

+ Y^{L,-{L,)md^,, (5.18) 

i 

where if is a Hamiltonian (of the open system) and Lj are Lindblad operators, as before, 
{Li) = {ip\Li\ip). The complex Wiener processes d^i satisfy 

M{d^i) = 0, M{diidQ = 0, M{di*dQ = 6ijdt, (5.19) 

where M denotes a mean over the ensemble. Quantum state diffusion reproduces the master 
equation in the mean: 
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p = M|V^)(V'|, (5.20) 



where IV') satisfy the quantum state diffusion equation ( [5.18| ), then it can be shown that p 
satisfies the master equation ( |5.17| ). 

In order to show the locahzation properties of the Langevin-Ito equation we now consider 
the simplest case which is assumed that bath temperature is zero (n = 0). In this case the 
master equation ( p.l|) reduces to 

-\- —{2apa) — a)ap — pa}a)^ (5-21) 
Then the corresponding Langevin-Ito equation is given by 

\d^) = -l.H\^)dt 



+|(2(a^)a - a^a - {a^){a))\i!)dt 



+Vl{a-{a)md^, (5.22) 

where is the complex Wiener process satisfying 

M{dO = 0, M{didi) = 0, M{dCdO = dt, (5.23) 

where M denotes a mean over probability distribution. 

The evolution of the quantum average of operators can be calculated by using the fol- 
lowing formula: 

d{G) = ^{[H,G]) 

~J2iLl[L^,G] + [G,L^]L,)dt 

i 

+ Y.i^{G\ U)di, + a(L,„ G)0, (5.24) 

i 

where 

a{A, B) = (A^B) - (A^) (B) (5.25) 
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Using Eq. (|5.24|) , it is straightforward to get the following equations: 



dt 



2 



1 + (a,) - + di 



d{ay) 



+ (o-^) - {a^f - i{a^){ay)\ dC, 

2 



(5.26) 



-^K)-yK) 



dt 



2 

V7 



(5.27) 



%/7 



2 

x/7 



{l + {a,)){{cr,) + i{ay))dC. 



(5.28) 



Moreover, we need to calculate the higher order moments. For any Hermitian operator A 



we have from (5.4) 



d{AAf = di{A^) - (Af) 

= d{A^)-2{A)d{A)-{d{A)f. 



(5.29) 



Then we easily obtain 



M—{Aa,)'' = 2u{a,){<yy) +^{a^f 

-|(Aa.)^-7(Aa.)2(a,) 
-|(a.)2-|((a,.)(a,))^ 
-2uj{a^){ay) +'y{cry)^ 



(5.30) 



7, 



-^(A^.) -7(Aa,)^(a,) 
-l{a^y-l{{a^){ay)y. 



(5.31) 



Now, we are in the position to consider the localization of solutions to Eq. ( p.22| ). Using 
master equation or quantum trajectories approach, it is very easy to see that the atom will 
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soon collapse into the lower state |1) and keeps there forever. Here we shall demonstrate 
that any solution to Langevin-Ito equation ( p.22| ) will localize at the lower state after a 
localization time. The collapse operator in this case is 



A = + icTy. 



(5.32) 



Then by using (W^) we get 



a{A, A) = {Aa^f + {Aayf + 2{a,). 



(5.33) 



Hence we have 



Mj^a{A, A) = 7(a.)2 + 7^)^ - ^(Aa.)^ 



-|7(Aa,)^ - 7(Aa.)^(a.) - 7(Aa,)2(a,) 
-7(a,)2-27((a.) + l)-7((a.)(a,))l 



(5.34) 



In order to prove that the left-hand side of Eq. ( |5.34| ) is non-positive, let us denote 



(Aa, 



(As 



1 + X, 



-1 + Z. 



(5.35) 
(5.36) 
(5.37) 



Substituting equations (|5.35|) , (|5.36|) and ( ^.37|) into Eq. ( ^.34|) we have 



Mj^a{A\A) = j 



-R^ - X -Y -2Z - ^{Y - Zf - ^(X - Zf 



(5.38) 



where R = {ax){<Jy). Note that 



X + Y + 2Z = a{A,A) > 0. 



(5.39) 



Then we show that 



M^a{A\A) < 

(JjL 



(5.40) 



and the equality holds if and only if 
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X = Y = Z = 0. (5.41) 

That is, the average in the left hand sides of equations ( |5.35| ), ( |5.3(j| ), and ( [5. 371 ) is taken 
over the ground state |1). This proves that the solution to Eq. ( |5.22| ) will localize at the 
ground state when the evolution time is larger than the localization time. 

Finally, let us estimate the localization rate of the quantum state evolution. Using Eq. 
( ^.38| ) and Eq. ( |5.39|) , we immediately obtain 

Mj^a{A,A)<-^iaiA,A)Y. (5.42) 



So the localization rate tiocaiization is 



^localization ~ — ■ (5.43) 

7 



Again we see that the above result is in agreement with the decoherence time scale (|2.19| ). 



VI. CONCLUDING REMARKS 

In this paper, based on the two-level system models, we have studied in detail the 
decoherent histories approaches. We present an explicit analysis of the degree of decoherence 
and its relation to the von Neumann entropy. We have demonstrated the localization in both 
quantum jump simulations and quantum state diffusion approaches. Here we conclude with 
a few remarks. 

Firstly, we have shown that the environment-induced decoherence, decoherent histories 
and the localization process in the stochastic Schrodinger equation approaches generally 
agree each other. We have shown that the set of histories consisting of the projection onto 
the basis in which the density operator tends to become diagonal give the better degree of 
decoherence. These results are in agreement with former studies on the quantum Brownian 



models [|I^ as well as the quantum optical models [Q. 

In this paper, we have shown that there are a number of sets of decoherent histories in 
this two levels model. Clearly, these decoherent histories are not equally important from 
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physical point of view. Among those, the most natural one is that which consist of the 
projections onto |1) and |2). We have proven that this set of histories give the best degree of 
decoherence. Note that the density matrix in the basis |1) and |2) will become diagonal after 
a typically short time. Moreover, we show that the solutions to the stochastic Schrodinger 
equation in the quantum jump simulation will localize at |1) or |2) after certain time which 
is basically same as the decoherence time. 

It is known that the solution generated by stochastic Schrodinger equation (|5.1| ) will 
randomly jump between the two levels |1) and |2). However, due to the influence of the bath, 
the average populations of the two levels will gradually become stable, and the localization 
process occurs. As we have shown, the set of histories consisting of the projection onto these 
two levels are perfectly decoherent. Here, we have seen that quantum jump simulation is 
entirely compatible with the history point of view. This is a very nice result. Similar results 



in the quantum state diffusion have been discussed before ||T7| , p!8 

In addition, we have found that the environment-induced decoherence, decoherent histo- 
ries and the localization process are more effective as the bath temperature increases. Phys- 
ically, this is an expected result as the bath at a higher temperature would have stronger 
influence on the system. As expected, the time scales concerning both decoherence and 
localization are basically same. 

Secondly, the approximate decoherence is of basic importance in practical physical pro- 
cess. By using this two-level system model we can clearly see what determines the degree 
of decoherence. For a given set of histories, the only adjustable parameters are the temper- 
ature of bath, the time-spacing interval and the initial state of the system. We also see a 
interesting relation between the degree of decoherence and the von Neumann entropy. This 
relation indicates that there is a tension between the predictability of state of system and 
the degree of decoherence. A similar tension has been discussed in a different situation based 
on quantum Brownian motion models |]^. 

Thirdly, it is important to notice that, as phenomenological theories, both quantum 
state diffusion and quantum jump simulation must be used under some conditions (e.g. 
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see |T2|). The comparison between different approaches therefore must be made in caution 
since the correspondence between them is by no means mathematically one-to-one corre- 
spondence. Rather, we emphasize that, underlying the open quantum system, the mutual 
influence between the system and its environment is the common theoretical base of all of 
those approaches and both decoherence and localization are nothing more than the different 
manifestations of a single entity. 

Finally, let us note that the coarse-graining in our two-level system is made by using 
projection operators on the system whilst ignoring the environment. It would be interesting 
to consider the general n-dimensional model in which the effect of a further coarse-graining 
on the degree of decoherence can be discussed. Work towards to this aspect is in progress. 

The environment-induced decoherence, decoherent histories as well as various stochastic 
Schrodinger equations have provided many important insights into the understanding of fun- 
damental problems in quantum theory. The investigation into the similarity and difference 
between the different approaches is of importance. The more thorough studies in this aspect 
would be useful. 
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APPENDIX A: PROOF OF THEOREM 

In this Appendix, we shall give a proof of Theorem in Section III. Since both M and N 
are positive definite matrices, therefore, one of them, say, can be decomposed as 
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= S^S, (Al) 

where S is an nxn matrix. After an arrangement, the right-hand side of Eq. ( [4 .31) becomes 

|Tr (^MPS^SQ) I = |Tr (^SQMPS^) \. (A2) 
Suppose orthonormal basis in n dimensional space V. Then 

I Tr (MPNQ) I = I ^ (^SQMPS^) x^l, (A3) 

m 

where is the transpose of Xm- Now, we set 

Vm = QS^Xm, (A4) 
= PS^Xra- (^5) 

Then the trace in Eq. ( [A^ ) may be rewritten as 

TT{MPNQ)=Y.ylMzm. (A6) 

m 

Since i/m, are orthogonal vectors and M is a positive definite matrix, then it is not difficult 
to arrive at the following inequality (see ||25|| ) 

I [fmMz^ I < e*^ [vl^My^ {z^Mzm) , ( A7) 

where = (A^,, - X^iJ/iX^,, + A^^jJ, A*^,, and X^n are the largest and the smallest 
eigenvalues of M, respectively. Combining (|A7D with Cauchy's inequality, 

(E«nM'<E«'E&n, (A8) 
m 

then (|A3|) becomes 



m m 



|Tr (MPNQ) I = I E^^m {SQMPS^) Xm\ 

m 

<J2\ylMzm\ 

m 

< e'\Y.ylMy^f'\T.^-lMz^f^\ (A9) 
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It is easy to identify that 

IV {MPNP) = Y: ylMy^, (AlO) 

m 

Tr {MQNQ) = ^ z^Mz^. (All) 

m 

This proves that 

|Tr (MPNQ) I < [Tr (MPNP)]^^^ [Tr (MQNQ)]^^^ . (A12) 

Since M and are in the completely symmetric position, so the similar result is true for 
e^. Then it completes the proof of the theorem. □ 
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